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ABSTRACT 

We look for a non-Gaussian signal in the WMAP 5-year temperature anisotropy maps 
by performing a needlet-based data analysis. We use the foreground-reduced maps 
obtained by the WMAP team through the optimal combination of the W, V and Q 
channels, and perform realistic non-Gaussian simulations in order to constrain the 
non- linear coupling parameter /nl- Wc apply a third-order estimator of the nccdlet 
coefficients skewness and compute the statistics of its distribution. We obtain —80 < 
/nl < 120 at 95% confidence level, which is consistent with a Gaussian distribution 
and comparable to previous constraints on the non-linear coupling. We then develop an 
estimator of /nl based on the same simulations and wc find consistent constraints on 
primordial non-Gaussianity. Wc finally compute the three point correlation function in 
nccdlet space: the constraints on /nl improve to —50 < /nl < 110 at 95% confidence 
level. 

Key words: cosmic microwave background - early universe - methods: data analysis. 



1 INTRODUCTION 



With the increasing amount of hl Kh-quality obse r vation s 
performed in th e las t decade ( Hinshaw et al.l (12008 



iReichardt et al 
2008t): 



(l2008l ): 



Prvke et al, 

l|2006h : " I Johnson etal 



^ Sieyers et al 
20081) 



(12007) 



Hinderks et al.l 

(|2007h ). statistical 



Wu et all 
2008t ): 



iMasi et al 

tests of the CMB temperature anisotropy pattern are 
getting more and more accurate. This has made possible to 
test one of the basic tenets of the standard cosmological sce- 
nario, i.e. that the primordial density perturbations follow 
a Gaussian distribution. This is a definit e prediction of the 
simp l est inflationary models (lGuthlll98ll : ISatolll98ll : iLindd 
ll982l : lAlbrecht fc Steinh ardt'> 19821): the detection of primor- 
dial deviations from Gaussianity would be a smoking gun 
for more complicated implementations of the inflationary 
mech a nism, such as those of mu l ti-fields (iLyth fc WaMg 
l2002l : iLinde fc Mukhano\j [200^: lAlabidi fc LythI l2006t ). 
ekpyrotic dMizuno et al.l 2008': 'Kh ourvl |2002|) or cyclic 
scena rios (jSteinhardt fc Turoh i2002l : iLehners fc Steinhardtl 
l2008l ). 



When dealing with the search for a non-Gaussian statis- 
tics in real data, two major issues have to be addressed. One 
has to do with the statistical tools used to analyse the data 
and detect deviations from Gaussianity: not only can these 
deviations be very subtle and elusive, but they could be gen- 
erated by processes that are not directly related to the pri- 
mordial density pertur bations — such as u nremoved astro- 
phys ical foregrounds (|Coorav et al ] l2008l : ISerra fc CooTwl 
l2008h or instrumental systematics. The other issue is the- 
oretical, and relates to the way the non-Gaussianity is pa- 
rameterised: while there is only one way to realize a Gaus- 
sian distribution, non-Gaussian statistics can be produced 
in countless ways. One then has to assume a non-Gaussian 
parameterisation which relates in some sensible way to an 
underlying early universe scenario. 

The latter issue is usually addressed by introducing 
a parameter /nl, which quantifies the amplitude of non- 
Gaussianity as a quadratic deviation with respect to the 
primordial Gaussian gravitational potential $l, i.e.: 



$(x) = $L(a:) + /nl ['I'l(x-) - {<i>Ux))] 



(1) 
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The major advantage of this parameterisation is that, re- 
gardless of the specific underlying early universe model, it 
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can represent the second-order approximation of any non 
linear deviation from Gaussianity. F or an excellent review 
on this topic see iBartolo et all (|2004l 'l. 

From the point of view of data analysis, a number 
of techniques have been proposed in the past few years 
to quantify the level of deviation from Gaussian statistics 
in the data. The most used one in harmonic space is the 
bispectrum (ILuoI Il9 94': 'Heavens' 'l998'; 'Spcrgcl fc Goldberd 
1 19991 : IKomatsu fc Spergeli i2001i : iCabcUa ct al.. ,200l ). The 
bispectrum is defined as the three-point correlation func- 
tion, and an estimate of /nl through the bispectrum re- 
quires the sum over all the triangle configurations. Since 
this is extremely time-consuming, regardless whether the 
computation is perfor med in harmonic space or pixel space, 
IKomatsu ei~^ (|2005l 'l have proposed a fast cubic estimator 
based on the Wiener filter matching, which reduces consid- 
erably the computati onal challenge. Th i s esti mator has been 
further developed bv ICreminelli et al.l l)2006h introducing a 
linear correction which allows the corre ct treatment of the 
anisotropic noise, and finally optimised (jYadav et al.l 200?!) 
and extended to polarisation measureme nts ( Yadav et al.l 
|2008| ). Recently, lYadav fc Wandeltl (|2008l ) applied the cu- 
bic estimator to the WMAP 3-year data, finding a detec- 
tion of a primordial non-Gaussian signal at more than 2.5 
sigma. An indication of a primordial non-Gaussian signal 
has been also found by the WMAP collaboration in the 
analysis of th e 5-year dataset, alth ough with a lower con- 
fidence level (|Komatsu et al.ll2008l ). An interesting discus- 
sion on optimal and sub-optimal estimators can b e found in 
ISmith fc Zaldarriag a (2006). Se for further details IYu fc Lul 
(|2008h and lBabichI (|2005l '). 



Con cerning the meth od s in p ixe l space, Ide Troia et al.l 



200 it): iReichardt et all (fiooi 



ICurto et al 
als to several 



(2008 



iHikage et all (|2006h : 
2003) applied Minkow ski f u nction - 
CMB datasets; fCa bcUa et al.1 (j2005l ) 
applied local curvatur e on WMAP 1-year data and 
iMonteserm et al] ()2006l ) developed scalar statistics using 
the Laplacian as a tool to test Gaussianity. Alterna- 
tive indicators base d on skewness and kurtosis have 
been proposed by iBernui fc Reboucad (|2008h . Tests 



(jVielva et al.ll2004: Mukheriee fc Wand 12004 


), and WMAP 


5-year data bv ICurto et al.l (20081) and 


McEwen et al.1 



Wavelets have b een a lso used in CMB studies 
2002h to i d entify anoma l ies in 



Martmez-Gonzalez et al 



Cruz et al. 2008a 


. 2007: Vielva et al. 20071: Wiaux et al. 


2ooe 


: Cruz et al.l 


20051: Vielva et al.l 20041: Pietrobon et al. 


2008 


), denoising \ 


Sanz et al.l 19991). point sources extraction 


(Cavonetal! 2000l: Gonzalez-Nuevo et al.ll2006l). Verv 


recentlv Vielva fc Sanz (2008h constrained primordial non- 



Gaussianity by means of N-point probability distribution 
functions. 

In this work, we constrain the primordial non- 
Gaussianity by applying for the first time a novel rendi- 
tion of spherical wavelets called needlets to the WMAP 
5- year CMB data. The n e edlet construction is discussed 
bvlNarcowich et all (|2006|) : iBaldi et al.1 ||200§): iBaldi et al.1 



(|2007^ : lMarinucci et al.1 (|2008l ): ICeller fc Marinuccil 
Needlets have a number of interesting properties which make 
them a promising tool for CMB data analysis: they live on 
the sphere, without relying on any tangent plane approxima- 



tion, are quasi-exponentially localised in pixel space and can 
be easily constructed from a filter function with a finite sup- 
port in harmonic space. Combined with the specific shape 
of this function, this guarantees a tiny level of correlation 
among needlets, which are only marginally affected by the 
pr esence of sky cut s . An e xhau stive discussion can be found 
in IPietrobon et al.1 (|2006l ) and iMarinucci et al.1 (|2008l ). Re- 
cently, the needlet formalism has been extend e d to t he the 
polarisation field, as discussed bv lGelle~ et al.1 (|2008l ). Here 
we try to test whether needlets can lead to interesting con- 
straints on the non-Gaussian amplitude, when compared to 
previous estimates of /nl- 

This paper is organised as follows: in Sec. [2] we review 
the needlets formalism; in Sec. |3] we describe the statistical 
analysis of WMAP 5-year data and the constraints on /nl ; 
we summarise our conclusions in Sec. |3] 



2 NEEDLETS FORMALISM 

We perform our analysis of the non-Gaussianity of WMAP 
5-year data by means of needlets. So far needlets have been 
successfully applied to the study of the CMB in the con- 
text of the detection of the Integrated Sachs- Wolfe effect 
fISW ) |Pietrobon et al.| [ 2006l ). the power spectrum estima- 
tion (|Cruz et al.ll2008b) and the study o f deviations from 
statistic al isotropv (IPietrobon et al.ll2008l ). We refer to the 
work bv IMarinucci et al.l ( 20081 ) for details on the construc- 
tion of a needlet frame and a detailed analysis of its sta- 
tistical properties. Here we remind that a set of needlets 
has one free parameter, B, which controls the width of the 
filter function in harmonic space. The filter function is non- 
vanishing only within a certain multipole range, thus allow- 
ing to keep trace of the angular information of the signal; 
the specific shape of the filter function also ensures a sharp 
localisation in pixel space. The correlation among needlets 
of the same set characterised by different j is much smaller 
than any other wavelet-like frame found in literature so far 
and can be easily described analytically. These features are 
particularly useful when looking for weak signals, like the 
non-Gaussianity one, which can be hidden by sky cuts and 
anisotropic noise. 

Formally, a needlet, ^lijki is a quadratic combination of 
spherical harmonics which looks like 



V'jfc(7) = 



(2) 



where 7 is a generic direction in the sky and ^jt refers to 
a cubature point on the sphere for the j resolution. The 
function b{£/B-') is the filter in ^ -space and A,fc is a cub ature 
weight number, which, foUowing lMarinucci et al.l(|2008h , has 
been set for practical purpose e qual to 1/Np, being Np the 
number of pixels in the Healpix (|G6rski et al.|[2005[ ) scheme 
for the resolution chosen. We show an example of the filter 
function for one set among those used in the analysis in 
Fig. [1] In the following we use the short notation bt for the 
filter function. For a given field defined on the sphere, T{-y), 
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b -filter 




1000 



Figure 1. Filter function in ^-spacc which the needlet construc- 
tion relies on. Set computed for B = 2. 



the needlet coefBcients are defined as 



(3) 



One of the most interesting properties of j3jk which makes 
needlets particularly suitable for handling CMB data is the 
following relation: 



jk 



4-K 



(4) 



which describes the square of needlet coefficients as 
an unbiased estimator of the angular pow er spectrum 
(|Pietrobon et al.ll2006l : iMarinucci et al.ll2008l ). This means 
that even if we group multipoles and each needlet peaks at 
a certain multipole range the total power is conserved: this 
property is peculiar of the needlets and it is not shared by 
other wavelet constructions. Note that this is strictly true 
only if mask are not applied to the maps, while small differ- 
ences are expected at low js if sky cuts are present. However 
when a rather broad and symmetric mask is applied, like in 
the case of Kq85 WMAP mask, the effect is very small. 

The signature of non-Gaussianity appears into the 
higher moments of a distribution, which are no longer com- 
pletely specified by the first moment (i.e. the mean value of 
the distribution) and the second moment (i.e. the standard 
deviation). For a Gaussian distribution, all odd moments are 
vanishing, while the even ones can be expressed in term of 
the first two only. We then look for a non-vanishing skew- 
ness of the distribution of the needlet coefficients, applying 
a cubic statistics. 



3 WMAP-5 NEEDLET ANALYSIS 

In the following we describe the statistical techniques and 
simulations used in order to constrain /nl. 

We started by producing simulations of non-Gaussian 
CMB maps with the WMAP-5 characteristics, for varying 
/nl. For each channel [Ql, Q2, VI, V2, Wl, W2, W3, W4], 
we used as input a realization of a simulated primordial 



non- Gaussian map (jLiguori et al.ll2007l ): these maps were 
convolved with the respective beam window functions for 
each channel and a random noise realization was added to 
each map adopting the nominal sensi tivities and number o f 
hits provided by the WMAP teanQ (|Hinshaw et al.1 120081 ). 
From th ese single-channel m aps we constructed an optimal 
map via IJarosik et all (|2007l ): 



r(7) =E^-'^W"'='>W 



(5) 



where 7 represents a direction on the sky (which, in prac- 
tice, is identified wit h a specific pixel in the Healpix scheme 
(|G6rski et al.l bOOSh ). and Wch = nh{-f)/ach/ J2ch'^':h- 
where nu is the number of observations of a given pixel 
and (Jch the nominal sensitivity of the channel, estimated 
by WMAP team. We finally applied the WMAP mask 
Kq85 and degraded the resulting map to the resolution of 
A'^ = 256. At the end of this procedure we were left with 
realistic Monte Carlo simulations of the CMB sky as seen 
from WMAP-5, containing difi'erent level of primordial non- 
Gaussianity parameterised by the value of /nl. 

Then, we extracted the needlet coefficients l3jk from the 
simulated maps for a given B. For each j resolution, the 
needlet coefficients can be visualised as a sky map, where 
k is the pixel number. We calculated the skewness of the 
reconstructed coefficients maps over the unmasked region, 
as: 



-^E 



(/3.fe'-(/3.fe'»' 



(6) 



where A'p denotes the number of pixels outside the mask 
and Oj is the variance of the needlet coefficients at the j 
resolution. This procedure allows us to build an empirical 
statistical distribution of the skewness as a function of /nl. 
Finally, we calculated the skewness from the real data of 
the foreground- reduced WMAP 5-year Q, V and W chan- 
nels data, using the same procedure applied to the simu- 
lated maps. The comparison of the real data skewness to 
the simulated distributions allowed us to set limits on the 
non-Gaussian signal in the data. 

A non- vanishing skewness represents a deviation from a 
Gaussian distribution and could give a fundamental handle 
on the physics responsible for inflation and the generation 
of primordial fluctuations. In general we expect the needlet 
coefficients to pick up signal at different angular scales as 
a function of both j and B, making different sets sensitive 
to non-Gaussianity in specific ways. This could be indeed 
a powerful tool when looking for imprint of specific models 
of non-Gaussianity. This feature is enhanced by the statis- 
tics itself we consider in our analysis. Since S oc we 
have an intrinsic modulation in the power of the cube of the 
needlet coefficients. Figure [2] shows this effect compared to 
the square of the hi function. 

For this reason it does make sense to compute the statis- 
tics defined in Eq.|6]for several sets of needlets. In particular 
we employed values of B from 1.8 to 4.5, choosing the step in 
order to span as homogeneously as possible the entire range 
of multipoles i" = 2 to £ = 500. The set of B we have consid- 
ered is [1.8, 1.9, 2.0, 2.15, 2.5, 3.0, 3.5, 4.0, 4.5]. We also tried 



http: / /lambda. gsfc.nasa.gov/product/map/dr3/m_products.cfm 
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bj sum 




Figure 2. Solid line sum of the dot-dashed Une sum of the 6|. 
While the former is equal to 1 for the entire range of multipole, 
the latter is not. 



Figure 5. The Ax^ of WMAP 5-year data as a function of /nl- 
/nl is estimated to be /nl = 20 ± 50 and /nl = 20 ± 100 at 1 cr 
and 2 cr level respectively. 



finer samplings of B, but no additional information resulted 
for the sampling considered. 

In Fig. |3] we report the skewness of the needlet coef- 
ficients, computed with the experimental set-up described 
above, for several values of _B as a function of the multipole 
£ where the resolution j peaks. The curves deviating from 
zero corresponds to the effect due to the primordial non- 
Gaussianity for positive (dashed lines) and negative (dot- 
ted lines) values of /nl, while the yellow and orange bands 
correspond to the la and 2a levels respectively. Diamonds 
represent the results of WMAP 5-year data. 



3.1 X analysis 

In order to estimate /nl we minimised the quantity: 

x'(/nl) = iX-" ~ {X{f^^))fC-\X^ - (X(/nl))). (7) 

Here X is a vector composed by the set of skewness of our 
set {B,j). The averages (X(/nl)) were calculated via Monte 
Carlo simulation over 100 primordial non-Gaussian maps. 
Formal ly, is dependent on /nl as well but it has been 



shown (|Spergel fc Goldberg|ll999i :[ Komatsu fc Spergelll200ll ) 
that for the mildly level of non-Gaussianity we expect this 
dependence is weak and can be estimated by Gaussian sim- 
ulations. We found that calculating C^^ from 10000 Monte 
Carlo simulations gives very stable estimate. 

The result is shown in Fig. O /nl is estimated to be 20 
with -30 < /nl < 70 and -80 < /nl < 120 at 1 and 2a re- 
spectively. These results show no significant deviation from 
the Gaussian hypothesis. This is not in contrast with the 
values found by Yadav fc Wandeltl (|2008l ). since we have per- 
formed our analysis on maps with the m aximum multipole 
corres ponding to £max = 500, whereas lYadav fc Wandeltl 
(|2008l ) clearly showed that their results crucially depend 
on the maximum multipole considered. The reduced of 
WMAP data is 1.53 with 85 degrees of freedom. 

As a further consistency check, we performed a 
goodness-of-fit test by calculating the quantile of the 
WMAP data from the non-Gaussian realizations with /nl = 
20. We found that 21% of the simulations had a larger x^ 
of the WMAP confirming that the specifications of our 
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Figure 6. Skewness for /nl = il, respectively dashed and 
dotted line, derived from non-Gaussian simulations. Theoretical 
curves for /nl = lOi 20, 30 (from bottom to top) are shown too. 



Monte Carlo simulations well describe the experimental set- 
ting of WMAP 5-year data. 



3.2 A skewness based /nl estimator 

Since the primordial non-Gaussianity is a second order ef- 
fect, it contributes linearly to the skewness through /nl. 
This means that it is possible to compute from the non- 
Gaussian simulations the skewness S(j) for /nl = 1 and 
use it as template to build a filter-matching estimator of 
the non-linear coupling parameter. Assuming that 5°''° = 



/nl Sj 



/nl- 



=1, where "th" means the average over the non- 



Gaussian simulations, we obtain 



(8) 



where we dropped the subscript /nl = 1. The theoretical 
skewness computed for /nl = 1 is shown in Fig. |B] We 
checked that the pipeline applied to simulated non-Gaussian 
CMB maps does not affect the linear relation: in particular 
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n = i.RO 




B = 2.00 



Figure 3. Skewness of needlet coefficients for different values of parameter B as a function of multipole i (logaritlimic function of tlie 
resolution j). Siiaded areas represents 1 and 2 sigma confidence levels calculated from 10000 Gaussian Monte Carlo simulations with 
beam, noise level, and mask of WMAP 5-year data. Dashed (dotted) lines correspond to the averages over 100 primordial non-Gaussian 
maps with /nl = 200,400 (-200,400). Diamond are the WMAP 5-year data. 



we verified that the average signal we obtain for a given /nl 
scales linearly with /nl itself, meaning for example that we 
can mimic the signal for /nl = ±400 by taking the double 
of that for /nl = ±200. 

The main contribution to the covariance matrix CoVjj/ 
comes from the Gaussian part of gravitational potential; this 
allows us to estimate the covariance from random Gaussian 
simulations. According to this assumption, we estimate the 
error bars on the primordial non-linear coupling parameter 
computing the standard deviation of the 10000 /nl esti- 
mates resulting from a fresh set of Gaussian simulations, 
via Eq. |8] We find /nl = 21 ± 54 at 1 sigma confidence 



level, which is fully consistent with what we found apply- 
ing the statistic. This corroborates the robustness of our 
procedure and confirms needlets as a suitable tool to study 
primordial non-Gaussianity. 

Our limits on primordial non-Gaussi a nity a re slightly 
larger than those achieved bv lCurto et all (j2008h bemg our 
1(7 confidence level (—30 < /nl < 70) slightly broader 
than —8 < /nl < —111 at 95% confidence level. However 
we were limited in our analysis to multipoles lower than 
imax = 500, while the strongest constraints on primordial 
non-Gaussianity make use of higher angular scales. 
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Figure 4. Skewness of needlet coefficients for tlie entire set of parameter B as a function of multipole £. Shaded areas represents 1 and 2 
a confidence levels calculated from 10000 Gaussian Monte Carlo simulations with beam, noise level, and marks of WMAP 5-ycar data. 
Dashed (dotted) lines correspond to the averages over 100 primordial non-Gaussian maps with /nl = 200, 400 (—200, 400). Diamond are 
the WMAP 5-year data. 



Our limits on /nl are quite promising for future ex- 
periments such as Plancfl where sensitivity and angular 
resolution will be enormously improved. 



3.3 Further analysis: binned bispectrum 

The analysis we performed in the previous sections, based on 
the skewness of the needlets coefficients Pjk, is mainly sensi- 
tive to the equilateral configurations, since it is proportional 
to the primordial bispectrum computed for £i ~ £2 — £3 
summed over all the multipoles up to £max = 500. This can 
be easily understood by direct inspection of Eq. |6l 



AT ^-^ 



« E E E bf^b^b'^a,,^,a, 

£imi £21^2 ^ams 

^£imi i^jk)y£2m2 fefc)^^3^3 fe^) 



27712 ^£3m3 



(9) 



2 http://www.rssd.esa.int/SA/PLANCK/docs/Bluebook-ESA- 
SCI(2005)l.pdf 



The sum over pixels approximates Gaunt 's integral 

k 
JS 

"(27i-H) (2^2 + 1) (2^3 + 1) 

£1 £2 £3 


£1 £2 £3 

mi m2 7713 



(10) 



and introducing the estimated bispectrum as the average 
over mi, m2 and 7713 



B 



£1 £2 



E 



m-l 771-2 771 



mi 1712 f^s 



we obtain the following relation: 



0), 0) , /(2^i + l) (2^2 + 1) (2^3 + 1) 



47r 



V , ^1 ^2 ^3 ,R 

X I )^'^'^'^ 



(11) 



Since the filter functions b'f ^ are computed for the same res- 
olution j, the biggest contribution to the bispectrum comes 
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Figure 7. The Ax of WMAP 5-year data as a function of /nl 
computed from the binned bispectrum. B is chosen 3.50 which 
had the highest signal-to-noise ratio among those chosen in the 
previous analysis, /nl is estimated to be /nl = 30 it 40 and 
/nl = 30 it 80 at 1 (T and 2 cr level respectively. 



from the equilateral configurations, plus a smaller effect of 
non-equilateral triangles when, for the same j, £i ^ £2 £3- 
A qualitative improvement in constraining the param- 
eter /nl can be achieved by adding in the estimator the 
effect of the squeezed configurations, considering the prod- 
uct of three Pjk with ji 7^ j2 7^ js. The skewness of needlet 
coefficients can be generalised into 



(12) 



'-2^3 

£1 £2 £3 



Ui ) ,U2 ) ,U3 ) . /(24 + l)(2^2 + l)(2£3 + l) 



47r 







can be seen as a binned bispectrum, a smooth and 



combined component of the angular bispectrum. 

We repeated the needlet analysis applying this new es- 
timator to the same set of WMAP 5-year data and simula- 
tions for the choice of the needlet parameter B = 3.5, which 
has the highest signal-to-noise ratio among the set chosen 
in the previous analysis. The minimization of the gives 
/nl = 30 ± 40, which is consistent with what we found ap- 
plying Sj. The for WMAP 5-year data is shown in Fig. [7] 

Recently an analysis based on a cu bic estimato r anal- 
ogou s to Sj-^j^jr^ has been performed bv lCurto et al.l (|2008l . 
I2OO9I ) who include the effect of squeezed co nfigurations in 
the Sp herical Mexican Hat Wavelet; and by iRudiord et aP 
(|2009l ) using a set of needlets characterised by a different B 
parameter obtaining /nl = 84 ± 40. While the difference in 
the value of /nl can be due to the higher number of mul- 
tipoles considered in the an alysis, ^max = 1000, and it is 
consistent with the result of lYadav fc Wandeltl (|2008l ). it is 
important that the estimated error bars are fully consistent 
with ours. 



4 CONCLUSIONS 

Primordial non-Gaussianity is becoming one of the keys 
to understand the physics of the early Universe. Several 
tests have been developed and applied to WMAP data to 
constrain the non- li near coupling parameter /nl. Recently, 
differ ent methods (jVadav fc Wandeltl bOOsI : ICurto et all 
|2008| ) brought to different constraints on /nl using simi- 
lar datasets, WMAP 3-year and WMAP 5-year respectively. 
The two approaches have been shown to have the same 
power in constraining primordial non-Gaussianity, while 
they obtained different best fit values for /nl. Whereas this 
might be due to the masks applied to the datasets, it cer- 
tainly underlines the complexity and difficulty of measuring 
/nl. The next generation of experiments will provide data 
with excellent angular resolution and signal-to-noise ratio 
which will be decisive to confirm or confute the measure- 
ments of /nl of the references above. In this respect, it will 
be even important to constrain /nl with different methods 
in order to get a more robust detection or to spot spurious 
presences of non-Gaussian signal. Moreover, integrated esti- 
mators, not based on Wiener filters, are differently sensitive 
to the non-linear coupling and can be useful to address ex- 
otic non-Gaussian models which predict high values of /nl 
and whose bispectrum evaluation, for instance, may require 
prohibitive computational time due to the convolution in 
the bispectrum formula. 

In this work, we constrained the primordial non- 
Gaussianity parameter /nl by developing the needlets for- 
malism and applying it to the WMAP 5-year CMB data. 
We estimated /nl to be 20 with —30 < /nl < 70 and 
— 80 < /nl < 120 at 1 and 2 sigma respectively, then consis- 
tent with the Gaussian hypothesis. We performed two dif- 
ferent analyses, the statistics and an estimator based on 
the skewness of the primordial non-Gaussian sky, finding 
an excellent agreement between the two results. Needlets 
have been proven to be a well understood tool for CMB 
data analysis, sensitive to the primordial non-Gaussianity. 
Since the skewness of the needlets coefficients is mainly sen- 
sitive to the equilateral triangle configurations, we improved 
our estimator computing the three point correlation func- 
tion in needlet space which indeed recovers the signal due 
to squeezed triangle configurations. We obtain /nl = 30±40 
at 68% confidence level, consistent with the previous analy- 
sis. Ouj_constraints_a^ slightly broader than those achieved 
bv lCurto et al. (l2008tl and not in contrast with the values 
found bv lYadav fc Wandel^ (|2008 ) since we were limited by 
a smaller range of multipoles, whereas the tighter constraints 
on /nl crucially depend on the maximum multipole consid- 
ered. 

Our limits on /nl are quite promising for future exper- 
iments like Planck, whose sensitivity and angular resolution 
will be enormously improved. 
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